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Abstract 

We perform the leading logarithmic resummation of nonglobal logarithms for the single-hemisphere jet 
mass distribution in e + e _ annihilation including the tinite-A> corrections. The result is compared with 
the previous all-order result in the large-iV c limit as well as fixed-order perturbative calculations. 


1 Introduction 

The hemisphere jet mass distribution is an event shape variable in e + e~-annihilation defined 
as the distribution of invariant mass in a single hemisphere whose axis coincides with the thrust 
axis. As is usually the case with all event shapes, it receives logarithmically enhanced pertur¬ 
bative corrections when the shape variable becomes small. However, unlike other event shapes 
for which systematic resummation methods are available (see (ll2l and references therein), the 
resummation of logarithms for the hemisphere jet mass distribution has turned out to be thorny 
and so far remained unsatisfactory even at the leading logarithmic level. This is due to the pres¬ 
ence of the so-called nonglobal logarithms (3] which arise from the energy-ordered radiation of 
soft gluons in a restricted region of phase space. 

The difficulty of resumming nonglobal logarithms stems from the fact that one has to keep 
track of the distribution of an arbitrary number of secondary soft gluons emitted at large angle. 
(For a recent review, see flU.) The original work by Dasgupta and Salam 0 employed a Monte 
Carlo algorithm, valid to leading logarithmic accuracy and in the large- N c approximation, to 
actually generate soft gluon cascades on a computer. Fater, Banff, Marchesini and Smye (BMS) 
(51 reduced the problem, still at large-iV c , to solving a nonlinear integro-differential equation. 
This latter approach paved the way for the inclusion of the finite- N c corrections in the resumma¬ 
tion (61 which has been recently put on a firmer ground (7J, and the first quantitative finite- N c 
result can be found in (8J. However, so far only one particular observable (‘interjet energy flow’) 
has been computed and the full impact of the finite- N c resummation is yet to be uncovered. 
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In this work, we apply the method developed in I® to the hemisphere jet mass distribu¬ 
tion and numerically carry out the resummation of nonglobal logarithms at fmite-iV c , thereby 
achieving the full leading logarithmic accuracy for this observable!!] In the next section we de¬ 
fine the observable and introduce the BMS equation which resums the nonglobal logarithms in 
the large-A^ c limit. In Section 3, we discuss the resummation strategy at finite- N c . It turns out 
that a naive application of the previous method is plagued by large numerical errors, and we 
shall propose a refined method which cures this problem. In Section 4, we present the numeri¬ 
cal result and compare it with the previous all-order result at large-iV c [|3]| as well as the recent 
fixed-order calculations fl9l l~0ll . 


2 Hemisphere jet mass distribution 


Consider a two-jet event in e + e -annihilation with the center-of-mass energy 0. Without loss 
of generality, we assume that the quark jet is right-moving with momentum = ^(1, 0, 0,1) = 
®(1, n R ) and the antiquark jet is left-moving with = ^(1, 0, 0, —1) = y(l, n L ). Suppose 
soft gluons with momentum fcf = cuj(l, n*) (n = (sin 6 cos </>, sin 6 sin (j>, cos 9)) are emitted in 
the right hemisphere 0 < < |. The invariant mass in the right hemisphere is 

M R = (p + E k i] = Q^ii 1 - COS0i) . (1) 

We shall be interested in the probability 

= ( 2 ) 

a Jo dp' 

that the rescaled invariant mass 
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(3) 


is less than some value p < 1. When p 1, one has to resum large logarithms In" 1/p in the 
perturbative calculation of Plr{p)- As observed in 0, this resummation consists of two parts. 
One is the Sudakov double logarithms (a s In 2 l/p) n which can be resummed via exponentiation. 
The other is the nonglobal logarithms (a s In 1 / p) n which arise from the fact that measurement is 
done only in a part of phase space (i.e., in the right hemisphere). The latter resummation affects 
various single-hemisphere event shapes in e + e“ annihilation and DIS j3). It is also relevant to 

1 Recently, the resummation of nonglobal logarithms to next-to-leading logarithmic order has been dis¬ 
cussed H7ll ill . In particular, Ref. Q explicitly derived the full NLL evolution kernel at finite N c . See, 
also, an earlier suggestion Ifl2l that the NLL resummation of nonglobal logarithms should be related to 
the NLL BFKL resummation via a conformal transformation. 
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the so-called soft function in the dijet mass distribution d 2 a/dM R dM L in the asymmetric limit 

m l » m r mim. 

So far, the resummation of nonglobal logarithms for P LR has been carried out only in the 
large- A c limit 03, or to finite orders of perturbation theory at large- N c 01 and finite- A c dna. 
As stated in the introduction, we shall perform the all-order leading logarithmic resummation 
at finite-along the lines of Il6l8il . To explain our approach, it is best to start with the BMS 
equation which resums both the Sudakov and nonglobal logarithms in the large-A c limit |(5| . 
Adapted to the hemisphere jet mass distribution 01, the equation reads 

drPap = N C J [Q L {l)P al P lP - Pap) , (4) 


where 


i A / \ 1 - cos 9 a0 

AI <2/3(7) /., p W-, p \ ’ 

(1 — COSt^JfT — COS U 1 p) 

is the soft gluon emission kernel and we defined 


(5) 


/ - -L-L-L . / 

7r p 

®l/r{i) is the ‘step function’ which restricts the angular integral dRL, = d cos 6~,d(t>~, to the 
left/right hemisphere. [Below we also use a shorthand notation f L / R d£l to represent this.] In 
©, P a/ 3 = P{^l a i O/^) is the generalization of Plr defined for arbitrary pairs of solid angle 
directions. 

Taken at its face value, Eq. © is ill-defined. When a or (3 is in the right hemisphere, the 
integral in the second term on the right-hand-side (the virtual term) is divergent, and this is 
precisely the situation (a/3) = [LR) we are eventually interested in. Physically, this collinear di¬ 
vergence should be cut off by the kinematical effect, yielding the Sudakov factor e ~ ln2 [ /f\ 
However, since the Sudakov factor is well understood anyway, one can leave it out of consider¬ 
ation by defining 


Pa /3 = exp y— 2 CpT j ^ ^.A/f 073(7) j 9 ap ■ (V) 

(C F = 1 and 2Cj? ^ N c in the large-A c limit.) Unlike ©, the equation satisfied by g a/3 is 

well-defined and amenable to analytical and numerical approaches. In particular, Ref. J9l ana¬ 
lytically calculated g R R to five loops using the hidden SL(2,M) symmetry of the BMS equation 

m. 


3 








3 Resummation at finite N c 


We now turn to the physically relevant case N c = 3. Temporarily forgetting about the issue of 
the collinear divergence, we recapitulate the resummation strategy developed in 111616181 . First 
we make the formal identification 

Pa/3 TT tr U a Up , (8) 

where U a is the Wilson line in the fundamental representation of SU(JV C ) from the origin to 
infinity in the fl a direction. The product in © represents the propagation of the qq jets (‘dipole’) 
in the eikonal approximation. As r is increased, more and more soft gluons are emitted from 
the dipole and also from the secondary gluons. This can be simulated as a stochastic process in 
which the Wilson lines receive random kicks in the color SU(JV C ) space, and is described by the 
following Langevin equation in discretized ‘time’ r 00 

U q (t + s) = e iS ° > e <A "U a (T)e iB °e iS " > , (9) 


where 

s» = 

A a = 

Bn, = 



(* = 1 , 2 ) 


and £0), are the Gaussian noises 


GO) 


( 11 ) 


( 12 ) 


(eS l M e$V)> = s“s r , T ,s(n y - Siy )s ab s u . (13) 

This is equivalent to the following ‘Fokker-Planck’ equation to be compared with © 

d T (p al> ) t =N c j - (P a „)() - ^e«(7)<fWs}. ( 14 ) 

2 We write the evolution ([9]) in a slightly different, but equivalent form compared to Ref. I®. It should 
be understood that various exponentials are meaningful only to 0(e) ®, although in practice we keep 
all orders in y/e in order to preserve the unitarity of U a . 
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where (... denotes averaging over the noises. In principle, Plr ( t ) at fmite-iV c can be evaluated 
by computing -^-trL^f/Jj for a given random walk trajectory with the initial condition U a (r = 
0) = 1, and then averaging over many trajectories. In this calculation, it suffices to define U a in 
the left hemisphere and at a single point a = R in the right hemisphere. 

However, this strategy does not apply straightforwardly to our presentproblem. ( Plr )§ 
quickly goes to zero due to the collinear divergence in the Sudakov factorjjj One may try to 
regularize the divergence by introducing a cutoff <5 and extract the finite part 


(, 9lr(t ))$ = lim exp 2 C f t 
5 —>-0 \ 




(15) 


Unfortunately, this does not work in practice because ( P[r )§ becomes very small and the ex¬ 
ponential factor becomes very large as 5 —» 0. It is difficult to numerically achieve the precise 
cancelation between the two factors. 

As a matter of fact, the same problem was already noticed in the original Monte Carlo sim¬ 
ulation at large-A^ c [O. There the authors subtracted the Sudakov contribution step-by-step, by 
modifying the emission probability as M. a p( 7 ) —* M. a p( 7 ) — Qr{i)M.lr{i)- Here we shall 
implement a similar subtraction directly in the evolution of U a . The origin of the collinear di- 
vergence can be traced to the factors e lSr > in ©. They give, after averaging over the noise, 


(e ba e ba e p e p )% 



f a ^) k 


^>7 a 


exp — 1 \ 


47t Jr 7 1 


(™/3 


n. 


np ■ n 0 


j-ac(i)k 
1 S7 a 


(16) 


which is indeed the Sudakov factor (0 generated in a single step. (fT6l) can be checked by using 
the identity [ 8 ] 


M a p( 7 ) = 2K, a p( r y) - £^( 7 ) - Kppffj) , K, a pin) = 


(n a - n 7 ) • (n 7 - rip) 


2(1 — n a ■ n 7 )(l — n 7 • np) 
In the special case (a/3) = ( LR ), we have that K/ LR (^) = 0 and 


.(17) 


Jr 17 A,1s(7) = /„ 17( _/C “ (7) “ Krr(7) ) = IT + L 




r Ati 1 — riR ■ n~ 


(18) 


We see that the singularity entirely comes from the second order term in the expansion of e 


iS 


(<) 


3 This problem was not encountered in |[ 8 ] because there a and /3 were always confined in the unob¬ 
served part of phase space. 
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It is thus tempting to remove the factor e a altogether and use a modified evolution equation 
U a (r + e) = e* Aa, I7 a (r)e* BQ . However, this also removes an essential part of the nonglobal 
logarithms. The reason is that the linear term in the expansion of e lSa = 1 + iS a (£) + • • • can 
give finite contributions when the Gaussian noise £ is contracted with that implicit in Uyf'UJ 
in (ITTI) . Physically, the factor U^t"UJ = (with Ua being the Wilson line in the adjoint 
representation) represents the emission of real gluons which is restricted to the left hemisphere. 
These gluons together with the original qq pair form a QCD antenna which coherently emits the 
softest gluon in the right hemisphere, thereby producing nonglobal logarithms. 

To make the last statement more concrete, we follow the evolution ([9]) analytically up to 
r = 2e (two steps) and collect the non-Sudakov contributions. We find 


{qlr{t))^ ~ 2C f N c t 2 


= -C F N r r 2 


dfL r dQ\ 


II 47t Jr 47t 
f dfL r dfl) 


(ICllW + ICrr (7)) (WA) + /C 7 fl (A) - /C 77 (A)) 


= —7T 


F iy c 

C f N c t 2 
12 


l 47T Jr 47t 


-Mlr{ 7 ) (A4l 7 (A) + A4 7j r(A) — Mlr(X)) 


(19) 


in agreement with the lowest order (two-loop) result Il3l9l FI In the first line of (fl9l) . the factors 
KJ Ll and JC jR come from the linear terms in e lSh ~ 1 + iS L and e lSlt ~ 1 + iS R , respectively. 
They both seem to be essential for obtaining the correct result. 

Importantly, however, the term /C /R (A) vanishes when integrating over the azimuthal angle 

0 A 


r2 TT 


d(pxJC.yR(X) = / d(f>\, 

Jo 


cos 6 *a — 1 — cos + cos 6 *a 7 
2(1 — cos6 | a)(1 — cos ^a 7 ) 


0 , 


( 20 ) 


where we used cos 6\ > 0 > cos # 7 . Moreover, by following the evolution © a few more steps, 
it is easy to convince oneself that the linear term iS R does not produce nonglobal logarithms to 
all orders because this term always reduces to factors like IC y ( n ) R (\) (after contracting with the 
n-th gluon emitted in the left hemisphere) and vanishes when integrating over 0\ in the right 
hemisphere. This observation brings in a major simplification in our resummation strategy. We 
can neglect the factors e R in d9J) for a = R and use the modified Langevin equation 


U R (r + e) = e iA «U R (r)e iBR . 


( 21 ) 


As for U a in the left hemisphere, we may continue to use the same evolution (©. Actually, 
we can make a slight improvement which speeds up the numerical simulation. The two inde¬ 
pendent noises defined in the right hemisphere always give identical contributions for the 
observable at hand. Therefore, we can eliminate one of them and use a modified equation 

1 At this order, we have to interpret 2c 2 = r(r — c) ~ r 1 to correct an error in iteratively solving a 
discretized differential equation. 
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( 22 ) 


C/.lr + £ ) = exp 


Note the factor of \/2. One can check that (|22h leads to the same equation (fl4l) for the product 
of two Wilson linesH] Using these Langevin equations, we finally compute the average 


= e r&l " 2 F (trC/ i (T)C/J(T)) s 


(23) 


The multiplicative factor in front subtracts the finite part of the Sudakov factor (fT8l) which is 
included in the evolution of Ur. 


4 Numerical results and discussions 


The numerical procedure is explained in Ref. 0 which we refer to for details. We discretize 
the solid angle 1 > cos 6 > — 1 and 2n > (f) > 0 into a lattice of 80 x 80 grid points and put 
a SU(3) matrix U a at each grid point on the left hemisphere cos 6 < 0. In addition, we define a 
SU(3) matrix Ur at a single point Ur in the right hemisphere. The Gaussian noise £ is randomly 
generated at all grid points and at each time step @ We then evolve U a and U R according to (l22l) 
and (12TT) . respectively, with e = 5 x 10~ 5 6 and the initial conditions U a = Ur = 1. As in the 
previous work (8], we observe large event-by-event fluctuations. In order to obtain a reasonably 
smooth curve, we typicallv have to average over C>(10 3 ) random walks. Fig. Q] shows the result 
from 2600 random walks]© In the same figure, we make comparisons with the following results 
in the literature: The blue line is a parameterization of the all-order Monte Carlo result in the 
large-limit by Dasgupta and Salam (DS) [0 


g DS (t) = exp 


(-CVAc 


7t 2 t 2 1 + (ar/2) 2 \ 
12 1 + (6r/2) c ) ’ 


(24) 


with a = 0.85A C , b = 0.86A C and c = 1.33. Here we set Cp ~ A c /2 = 1.5 which is what was 
actually used in 0. The black dashed line is a combination of the fixed-order analytical results 

5 We have checked numerically that (l22l) and © give equivalent results. The equivalence may not hold 
for more complicated observables. 

6 In |H1, the authors inadvertently assumed that the noise £ (at each time step) is independent of 0 at 
the degenerate points cos 0 = ±1. Fortunately, this was innocuous for the observable considered in 
|[8ll . However, for the present observable this causes a systematic error in the evaluation of the Sudakov 
integral (the first term of (flST) ) already at small-r because the integration region f cIQr includes the point 
cos 9 = 1 (which was not the case in |[8|l). In the present simulations, we fixed this problem by generating 
£ at different values of (j) independently also at cos 6 = ±1. 

' We also performed simulations with different discretization parameters (160 x 40 and 80 x 40 lattices, 
and s = HF 4 ) and found that the results arc consistent with each other. 
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by Schwartz and Zhu (SZ) to five-loop at large-A/ [|9j and Khelifa-Kerfa and Delenda (KD) to 
four-loop at finite N c ea 


, SZ-KD, 


7 r 


t) = 1-C f N c —t z + 


CfN^s 

6 


25 


13 


17 


V2 1 r 4 


120 


(-8C 2 F N*( 2 ( 3 + y CfNGs) t 5 , 


(25) 


where Cf = 4/3. Actually, the complete finite- AT C result at O(r') is not available, and the 
above formula is a well-motivated guess iflOl which reduces to the known result in the large- N c 
limit. Finally, the blue dash-dotted line is the following ‘resummed’ expression suggested by 
KD based on their four-loop result 


7r 2 r 2 


fn(resum)^ = exp _ CpN + 


C F N*( 3 T 3 7T 4 /25 


12 


6 


135 V 8 


( -CfK ± ) - 


,(26) 


with Cp = 4/3. Note, however, that at the moment it is not known whether the nonglobal loga¬ 
rithms actually exponentiate to all orders. 


We see that our result agrees very well with the most-advanced fixed-order result (1251) up to 
r < 0.5. Beyond this, the perturbative result quickly deviates and eventually blows up. It has 
been observed that higher loop contributions alternate in sign and converge rather poorly IIT71 . 
In addition, fixed-order results are numerically sensitive to the 1/Ay-suppressed corrections 
when r ~ 0(1). This can be partly remedied in the resummed formula (l26l ) . On the other hand, 
the all-order large- N c result (l24l) stays close to our curve up to r = 1. In fact, the difference can 
be partly accounted for by choosing C F = 4/3 in (l24l ) . which is what was actually suggested 
by DS as the likely functional form at finite-A/ and has been used for phenomenological pur¬ 
poses mum This is shown by the red line in Fig. |T] To correct the remaining difference, we 
independently determined a,b,c in (l24l) with C F = 4/3 and obtained 

a = 0.62 ±0.06, b = 0.06 ± 0.03, c = 0.37 ±0.04. (27) 


In conclusion, we have completed the resummation project for the single-hemisphere jet 
mass distribution initiated in [|3]] by including the finite-A/ corrections to all orders. We find 
that the finite- N c effect is numerically small, and this is consistent with the previous finding 
in [J8]|. However, it should be kept in mind that the observables calculated at finite-A^ c so far 
are defined in e + e~ annihilation where the two outgoing jets are represented by the product of 
two Wilson lines tr U a Ul. In hadron-hadron collisions, or in processes including hard gluons, 
one needs to consider the evolution of more complicated objects such as tr(U a UlUjUl) and 
tr(U a Up)tx(UjUl) (cf., Il20l ). The finitc-Av effects in the resummation of nonglobal logarithms 
for these observables have not been studied so far. 
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Fig. 1. ( 5 Li?(r)) at N c = 3 as a function of r obtained by averaging over 2600 random walks. The error 
bars indicate the standard error. Data points are plotted every 0.01/e = 200 random walk steps. Various 
curves are explained in the text. 
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